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Abstract 

Background: Relatively recent evidence indicates that ABCC2 transporters play a main role in the mode of action 
of Bacillus thuringiensis (Bt) Cry1A-type proteins. Mapping of major CrylA resistance genes has linked resistance to 
the ABCC2 locus in Heliothis virescens, Plutella xylostella, Trichoplusia ni and Bombyx mori, and mutations in this gene 
have been found in three of these Bt-resistant strains. 

Results: We have used a colony of Spodoptero exiguo (Xen-R) highly resistant to a Bt commercial bioinsecticide to 
identify regions in the 5. exiguo genome containing loci for major resistance genes by using bulk segregant analysis 
(BSA). Results reveal a region containing three genes from the ABCC family {ABBC1, ABBC2 and ABBC3) and a 
mutation in one of them {ABBC2) as responsible for the resistance of S. exiguo to the Bt commercial product and to 
its key Spodoptera-active ingredients, CrylCa. In contrast to all previously described mutations in ABCC2 genes that 
directly or indirectly affect the extracellular domains of the membrane protein, the ABCC2 mutation found in 
5. exiguo affects an intracellular domain involved in ATP binding. Functional analyses of ABBC2 and ABBC3 support 
the role of both proteins in the mode of action of Bt toxins in 5. exiguo. Partial silencing of these genes with dsRNA 
decreased the susceptibility of wild type larvae to both CrylAc and CrylCa. In addition, reduction of ABBC2 and 
ABBC3 expression negatively affected some fitness components and induced up-regulation of orylphohn and repots, 
genes that respond to Bt intoxication and that are found constitutively up-regulated in the Xen-R strain. 

Conclusions: The current results show the involvement of different members of the ABCC family in the mode of 
action of B. thuringiensis proteins and expand the role of the ABCC2 transporter in B. thuringiensis resistance beyond 
the Cryl A family of proteins to include CrylCa. 
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Background 

Genes coding for insecticidal proteins from the bacter- 
ium Bacillus thuringiensis (Bt) have been transferred to 
plants to protect them from the attack of insect pests 
[1,2] and the global impact of this strategy has been un- 
precedented in plant biotechnology. Unfortunately, and 
similarly with the history of chemical insecticides, in- 
sects have developed resistance to B. thuringiensis crystal 
proteins (Cry proteins) when exposed to either foliar B. 
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thuringiensis sprays or transgenic crops (Bt-crops) [3-5]. 
The long lasting benefits of B. thuringiensis insecticides 
and Bt-crops will mostly depend on the understanding 
of the mode of action of Cry proteins and of the bio- 
chemical and genetic bases of insect resistance. 

All B. thuringiensis crystal proteins (Cry proteins) which 
are toxic to target insects share general traits in their mode 
of action: solubilization in the midgut, activation by gut 
proteases, binding to membrane surface proteins (generally 
referred to as receptors') [6,7]. Susceptible insects have 
epitopes in their membrane receptors that are recognized 
by sequences of the Cry proteins. This is probably the main 
reason why Cry proteins are not harmful to non-target 
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insects [8]. The key role of membrane receptors has also 
become evident when studying the biochemical basis of in- 
sect resistance to Cry toxins. Most cases of high level re- 
sistance to CrylA proteins have been shown to be due to a 
decrease of toxin binding to midgut receptors [9]. 

Two approaches have been undertaken to characterize 
membrane proteins involved in the toxicity of Cry pro- 
teins. Formerly, the approach was biochemical, and it 
allowed for the characterization of receptors such as 
cadherin, aminopeptidase N (APN), and alkaline phos- 
phatase (ALP) [10], among others. The occurrence of 
such a variety of receptors has been explained by inte- 
grating them in a 'toxicity pathway that implies their 
participation in a sequential way [7]. In the most ac- 
cepted model, binding to cadherin is required for the 
Cry proteins to oligomerize and then bind to other re- 
ceptors, such as APN and ALP, which favor toxin inser- 
tion into the membrane. There is evidence supporting 
the participation of cadherin and APN in the toxicity of 
Cry proteins [11,12]. Furthermore, mutations conferring 
resistance to proteins from the CrylA family have been 
found in the cadherin gene from Heliothis virescens [13], 
Helicoverpa armigera [14,15] and in Pectinophora gossy- 
piella [16]. Similarly, mutations or reduced expression of 
apn genes have been found in resistant insects from 
Spodoptera exigua [17], Diatraea saccharallis [18] and 
Trichoplusia ni [19]. 

Another membrane protein, identified as an ABCC2 
transporter, was found related to the mode of action of 
CrylA-type proteins using a genetic approach to identify 
major genes involved in resistance in H. virescens [20]. 
Later on, mutations in the ABCC2 gene were found in 
CrylA resistant colonies from three other insect species: 
Plutella xylostella [21], T. ni [21], and Bombyx mori [22]. 
Further evidence for the involvement of the ABCC2 
transporter in CrylA toxicity came from functional stud- 
ies with transfected cultured cells [23]. Although muta- 
tions in either the cadherin or in the ABCC2 gene alone 
may suffice to render high levels of resistance, it seems 
clear that the occurrence of both functional genes is ne- 
cessary for the full toxicity of CrylA proteins [20,23] . Al- 
though the binding of CrylA proteins to the ABCC2 
transporter has never been directly shown, it is thought 
that the ABCC2 transporter transiently binds the Cry 
oligomer and mediates its insertion from the APN or 
ALP receptors to the membrane [24]. 

In a previous work we have reported an S. exigua 
strain (Xen-R) with high levels (> 1,000-fold) of resist- 
ance to a B. thuringiensis -based pesticide, Xentari™ [25]. 
Expression of about 600 midgut genes was analyzed by 
DNA-macroarray in order to find differences in midgut 
gene expression between susceptible and resistant in- 
sects. Among the differentially expressed genes, repat 
and arylphorin were identified and their increased 



expression was correlated with B. thuringiensis resist- 
ance. We also found overlap among genes that were 
constitutively over-expressed in resistant insects with 
genes that were up-regulated in susceptible insects after 
exposure to Xentari™, suggesting a permanent activation 
of the response to Xentari™ in resistant insects. However, 
we could not determine the exact resistant gene or genes 
and the molecular mechanism responsible of the high 
levels of resistance to B. thuringiensis in these insects. 

In the present paper we used bulk segregant analysis 
(BSA) based on high-throughput sequencing methodolo- 
gies (also known as next generation sequencing' or NGS) 
[26,27] to identify regions in the S. exigua genome con- 
taining loci for major resistance genes in Xen-R strain 
[25]. The results reveal two ABCC genes (ABBC2 and 
ABBC3) involved in the mode of action of B. thuringiensis 
toxins and a mutation in one of them (ABBC2) as respon- 
sible for the resistance of S. exigua to a B. thuringiensis 
commercial product and to its main active ingredients, 
CrylCa and CrylA toxins. The current results expand the 
role of the ABCC2 transporter in B. thuringiensis resist- 
ance beyond the CrylA family of proteins to include 
CrylCa. 

Results 

Construction of backcross families for resistance mapping 

Mapping of resistance was approached by mixing the re- 
sistant and susceptible genomes followed by selection for 
resistance after one round of backcrosses. At the concen- 
tration of the Xentari™ product used, Xen-R showed in- 
complete resistance because 65% of the larvae did not 
survive the treatment. The same treatment killed 100% of 
FRA and Fl larvae, indicating that resistance was com- 
pletely recessive at the diagnostic concentration. After 
crossing Fl and Xen-R adults and selection with the dis- 
criminant concentration of the Xentari™ product, three 
families resulting from independent backcrosses were ran- 
domly selected for the BSA. Mortality in the treated sam- 
ples was 87% for BC1, 89% for BC3 and 84% for BC4. 
Considering that the baseline mortality for Xen-R was 
65%, the mortality in the three backcross (BC) families is 
in good agreement with resistance being due to one major 
locus. Survivors from the three BC families were pooled 
(Sel-BC sample) and used for the BSA. 

New generation sequencing, assembly and mapping 

Parental samples were differentially tagged and sequenced 
in a single Hi-Seq2000 paired-end (PE) line and Sel-BC in 
a second line. A summary of the sequencing and assem- 
bling process is reported in Table 1. These reads (GenBank 
SRA accession: SRS514057, SRS514097, and SRS514098), 
along with the 445,524 sequences from the 454-sequencing 
and the previously published 6,212 Sanger sequences [28], 
were assembled using Trinity, yielding a total of 138,989 
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Table 1 NGS results for the different samples and transcriptome statistics 



Sample 


Raw reads 


Cleaned reads 


% mapped sequences 


Xen-R 


179,215,722 


168,518,208 


88.57 


FRA 


200,970,866 


187,779,324 


88.19 


BC-sel 


388,370,500 


360,354,916 


87.27 


Transcriptome 


Non-filtered 


Filtered 




Unigenes 


138,989 


96,676 




Gene clusters 


85,81 1 


54,478 




N50 


1,781 


2,016 




Size (Mb) 


119 


109 




Drosophila orthologs 




10,421 




Silkworm orthologs 




15,116 




NGS, next generation sequencing. 



unigenes and 85,811 gene clusters with an average length 
of 881 bp, with a total length 119 Mb. These unigenes were 
filtered to remove small, non-highly expressed anonymous 
unigenes. Finally, 96,676 unigenes and 54,478 gene clusters 
were kept, with an average size of 1,124.25 and a total 
length of 109 Mb (GenBank BioProject 219058). These 
unigenes were compared using bi-directional Blast search 
with Bombyx and Drosophila proteins. For 10,421 and 



15,116 unigenes, orthologs were identified in Drosophila 
and silkworm, respectively (Table 1, and Additional file 1). 
All reads were mapped against a filtered reference tran- 
scriptome built with the longest unigene of each gene clus- 
ter. The proportion of the total sequences which could be 
properly mapped was 87.27%, 88.19% and 88.57% from 
Sel-BC, FRA and Xen-R, respectively. The average cover- 
age was 763.9x. 




&*P 9TJip 



Figure 1 Schematic representation of the BSA index for the SNPs in the S. exigua unigenes with orthologs into the silkworm genome. 

Each data point (from 0 to 1) is the average index for three consecutive orthologs calculated using overlapping rolling-windows (BSA index for 
each unigene individually is provided in Additional file 3). Average index >0.60 are painted in red. BSA, bulk segregant analysis; SNPs, 
single nucleotide polymorphism. 
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Bulk segregant analysis 

Using Varscan, 437,815 single nucleotide polymorphisms 
(SNPs) and 80,246 indels were identified [see Additional 
file 2]. To select unigenes with a biased allele distribu- 
tion in the Sel-BC sample, the following index was 
calculated: ratio of SNPs with the Xen-R allele fixed in 
Sel-BC (frequency higher than 0.95) against the total 
number of SNPs. This index ranged from 0 to 1. Unigenes 
with an index close to 1 had a bias towards the allele 
present in the Xen-R parental. This index was plotted 
along with the Bombyx genome (Figure 1) by using the 
position of the respective Bombyx orthologs. A total of 
978 unigenes [see Additional file 3]. With a silkworm 
ortholog were available. Two wide Bombyx genomic re- 
gions with a high bias towards the resistant allele were de- 
tected: a wider region in chromosome 15 and another 
smaller region in chromosome 22. The region in chromo- 
some 15 spanned from 64,363 to 3,068,277 of nscaf2888 
scaffold and contained homologs to 64 S. exigua unigenes. 
The region in chromosome 22 spread from 67,364 to 
1,716,412 of nscafl681 scaffold and contained homologs 
to 42 S. exigua unigenes. 

Detailed analysis of the chromosome 22 region did not 
reveal genes potentially involved in the mode of action 
of B. thuringiensis or other entomopathogens. Moreover, 
we could not identify any SNPs associated with Sel-BC 



or Xen-R insects having a major effect (for example, 
nonsense mutation) on the codified proteins when com- 
pared to the alleles in FRA insects. 

Chromosomal regions from several insect species hom- 
ologous to a region in chromosome 15 in B. mori have 
been previously associated with resistance to B. thuringien- 
sis CrylA toxins [20-22]. Such a B. mori region contains 
several members of the ABCC family. Detailed analysis of 
unigenes from S. exigua in that region identified three 
ABCC genes, which we named ABCC1, ABCC2 and 
ABCC3 according to a phylogenetic analysis with repre- 
sentative members of ABC transporter genes from other 
Lepidoptera and Coleoptera (Figure 2). The S. exigua 
ABCC2 and ABBC3 genes codify for a 12-transmembrane 
domain protein [see Additional file 4: Figure SI]. Although 
our transcriptome did not contain the complete ABCC1 
unigene sequence, comparison with other ABCC1 genes 
suggests that it codifies for a protein with a topology simi- 
lar to that of other ABCC1 transporters. Allelic compari- 
son of unigenes in that region revealed a major mutation 
in the ABCC2 gene: a deletion of 246 nucleotides (82 
amino acids) in the unigene from the Xen-R and Sel-BC 
samples, but not from the FRA sample. Despite the high 
coverage for that gene in the Sel-BC sample, no single read 
was observed mapping into the deleted region. The deleted 
region included a fragment of the predicted ATP binding 




ABCC2 



Figure 2 Phylogenetic analysis and ortholog asignment of S. exigua ABCC transporters. Neighbor-joining consensus tree of representative 

members of ABC transporters from Lepidoptera and Tribolium costoneum (Coleoptera). Nomenclature for each member is based on the name of 

the species (Brno: Bombyx mori; Se: Spodoptera exigua; Dpe: Danaus plexippus; Bma: Bombyx mandarina; Hvi: Heliothis virescens; Hsu: Heliothis 

subflexo; Mse: manduca sexta; Tea: Tribolium costoneum) followed by their accession number from the National Center for Biotechnology 

Information (NCBI), Silkdb and Manduca Base (Agripest) databases. Sequences were selected based on their homology to the ABCC1-3 

transporters from the linkage region in B. mori. Some homologs in T. cotoneum were selected as a kind of outgroup. For the sake of clarity 

bootstrap values are only reported for the branches identifying the main ABCC subgroups. NCBI, National Center for Biotechnology Information. 
\ J 
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domain 2 in the hypothetical protein (Figure 3C). Further 
analysis of the genomic region revealed a deletion of about 
500 nucleotides extending over two exons (Figure 3). 

Based on these results, we decided to focus on the 
possible function of the ABCC2 and ABCC3 proteins, 
especially on their role in the susceptibility to B. thurin- 
giensis toxins and the effect of the ABCC2 mutation on 
the binding of CrylCa toxin to the insect midgut. 

Expression profiles of SeABCC2 and SeABCC3 genes in 
S. exigua 

Gene expression of SeABCC2 and SeABCC3 was ana- 
lyzed by reverse transcriptase-polymerase chain reac- 
tion (RT-PCR) in different developmental stages of S. 
exigua and in different tissues of fifth instar larvae. 
SeABCC2 and SeABCC3 were expressed in all larval 
instars and in the adult stage although at early stages 
SeABCC2 was apparently more strongly expressed than 
SeABBC3 (Figure 4A). Regarding tissue distribution 
(Figure 4B), SeABCC2 was highly expressed in the gut and 
hemocytes; in addition, the sample from the fat body gave 
a very faint amplification band. SeABCC3 was also highly 
expressed in the gut and hemocytes and faint amplifica- 
tion bands were obtained in samples from the fat body 
and nerve, but not from the salivary glands. 



Effect of SeABCC2 and SeABCC3 suppression on the 
susceptibility of S. exigua to Cry! Ac and CrylCa 

To study the role of these proteins in the susceptibility 
of S. exigua to B. thuringiensis toxins, SeABCC genes 
were silenced by RNA interference (RNAi) and the sup- 
pressed larvae were tested for their susceptibility to two 
different Cryl toxins. Prior to performing bioassays with 
double-stranded RNA (dsRNA) -treated larvae, the sus- 
ceptibility of non-treated S. exigua third instar larvae to 
Cryl Ac and CrylCa pro toxins was determined under 
our experimental conditions. As expected, CrylCa was 
significantly more toxic than CrylAc, with LC 50 values 
of 0.24 and 4.8 (ig/ml, respectively. The toxin concentra- 
tions that yielded mortality between 90% and 100% were 
10 and 1 (ig/cm 2 for CrylAc and CrylCa, respectively 
(Figure 4C). These toxin concentrations were chosen for 
subsequent experiments with dsRNA. 

Expression of SeABCC2 and SeABCC3 started to de- 
crease after 12 hours of dsDNA ingestion, but the sup- 
pression was more drastic after 48 hours and lasted at 
least up to 72 hours in both cases (Figure 4D-E). Cross- 
suppression between the two genes was discarded by es- 
timation of SeABCC2 expression in ds«SeABCC3-treated 
insects and vice versa (Figure 4D-E). When S. exigua lar- 
vae were treated with either dsSeABCC2 or dsSeABCC3 
and then exposed to 10 [ig/cm 2 of CrylAc or 1 (ig/cm 2 





Figure 3 Localization of the deletion in resistant insects into the ABCC2 gene genomic region (A) schematically represented at the 

genomic region (B) and on the hypothetical protein (C). Sequences were obtained by PCR amplification of the reported region using 

genomic DNA from the Xen-R and FRA strains. 'E' labels are indicative of the predicted exons in the reported region. The dashed line indicates 

the deletion in the Xen-R insects. PCR, polymerase chain reaction. 
\ J 
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Figure 4 Expression patterns of SeABCC2 or SeABCC3 gene in different developmental stages (A) and in different tissues (B) of the fifth 
instar larvae in S. exigua. 'E', 'L1 - L5' ; V, and 'A' represent egg, larval instars, pupa, and adult, respectively. 'FB', 'HC, 'NV, and 'SG' represent fat 
body, hemocytes, nerve, and salivary gland, respectively. Expression of (3-actin confirms the integrity of cDNA preparation. Toxicity comparison of 
CrylAc and CryICa protoxins against third instar larvae of S. exigua (C). Each bioassay was conducted with ten larvae per replicate and three 
replicates per concentration. Mortality was scored at five days post-treatment. Suppression of SeABCC2 or SeABCC3 gene expression using a 
specific double-stranded RNA {6sSeABCC2 or 6sSeABCC3) in the whole body of S. exigua larvae at different periods (D and E, respectively). SeABCC2 
or SeABCC3 expression was analyzed by qRT-PCR. Each treatment was independently replicated three times. Different letters above standard 
deviation bars indicate significant differences among means at Type I error = 0.05 (LSD test). The effect of suppression of SeABCC2 or SeABCC3 
gene expression on susceptibility of the third instar S. exigua to CrylAc or CryICa protoxins (F). A viral gene (ORF302)-specific dsRNA (dsCON) 
was used for a dsRNA control. At 24 hours after dsRNA treatment, the larvae were exposed to CrylAc (10 ug/cm 2 ) or CryICa (1 ug/cm 2 ) protoxins 
treated by the leaf-dipping method. Mortality was assessed at five days after Cry toxin treatment. Each treatment was replicated three times. Each 
replication used 10 larvae. Different letters above standard deviation bars indicate significant difference among means at Type I error = 0.05 (LSD 
test). dsRNA, double-stranded RNA; LSD, least significant difference; qRT-PCR, quantitative reverse transcriptase-polymerase chain reaction. 



of CryICa, they showed a significant decrease in suscep- 
tibility to these two toxins (Figure 4F). In particular, si- 
lencing of SeABCC2 highly decreased CryICa toxicity to 
S. exigua compared with silencing of SeABCC3. Under 
these conditions, control larvae suffered about 92% and 
95% larval mortality to CrylAc or CryICa, respectively. 



Binding of CryICa to brush border membrane vesicles 
from resistant and susceptible larvae 

Since mutations in the ABCC2 transporters have been 
proposed to confer resistance to insects from different 
species that lacked binding of one or more CrylA pro- 
teins [20,21] binding of 125 I-labeled CryICa was tested 
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in susceptible and resistant insects to see whether bind- 
ing was altered in the latter (Figure 5). Competition 
binding curves indicated that resistant insects had a high 
contribution of non-specific binding (the binding that 
cannot be displaced with high concentrations of com- 
petitor) to the total binding (Figure 5A). The equilibrium 
binding parameters did not differ significantly between 
the two types of insects and the only significant differ- 
ences were found in the concentration of binding sites 
(Table 2). Specific binding of Cry toxins to their mem- 
brane receptors consists of a reversible binding compo- 
nent and an irreversible binding component associated 
with the toxin insertion into the membrane [29]. To 
search for additional possible differences in binding, an 
experiment was designed to dissect the total specific 
binding into its reversible and irreversible components. 
The results showed a significant decrease in the irrevers- 
ible component of the specific binding in resistant in- 
sects as compared with susceptible insects (Figure 5B). 
However, the reversible component was similar in both 
cases. 

Fitness-cost associated with the silencing of SeABCC2 or 
SeABCC3 expression 

The effect of silencing of either SeABCC2 or SeABCC3 
on the fitness of the insects was also measured. Reduc- 
tion in the expression level of both genes had a statisti- 
cally significant influence on both pupation and adult 
emergence. Feeding dsRNA specific to either SeABCC2 
or SeABCC3 reduced the percentage of pupation to 70% 
and 81%, respectively (Figure 6A) and the percentage of 
adult emergence to 72% and 68%, respectively (Figure 6D). 
Treatment of dsSeABCC2 also had a significant effect on 
the pupal weight, although that of dsSeABCC3 did not 
(measured at the second day of pupal stage) (Figure 6B). 
In contrast, none of the treatments had a significant effect 



Table 2 Equilibrium dissociation constant (K d ) and 
concentration of binding sites (R t ) of CryICa binding to 
BBMV from Xentari™ susceptible (FRA) and resistant 
(Xen-R) 5. exigua larvae 



Colony 




Mean ± SEM° 




K d (nM) 


R t (pmol/mg) 


FRA 


14 + 2 


5.9 ±0.6 


Xen-R 


28+13 


2.85 ±1.13 



a From two independent replicates. BBMV, brush border membrane vesicles; 
SEM, standard error of the mean. 



on the duration of the pupal stage, which in all cases took 
five to six days from the first day of pupation till adult 
emergence (Figure 6C). 

Overexpression of arylphorin and repatS genes induced by 
SeABCC2 and SeABCC3 suppression 

Since both arylphorin and repatS have been reported to 
be overexpressed in the Xen-R insects [25], it was of 
interest to determine whether the silencing of SeABCC2 
or SeABCC3 could have an effect on the expression of 
these genes. The results show, in agreement with previ- 
ous reports, that when either the SeABCC2 or SeABCC3 
gene is silenced, both arylphorin and repatS genes in- 
creased their expression as compared with that of the 
controls [see Additional file 5: Figure S2]. 

Discussion 

Forward genetics has been useful in shedding light on 
new resistance mechanisms to CrylA proteins in insects. 
Relatively recently, the role of the ABCC2 transporter 
has been made evident in the mode of action of CrylA 
proteins through linking resistance to mutations in this 
membrane protein in four different insect species [24]. 
We have successfully adapted standard BSA procedure 
to our particular conditions: 1) since we could not easily 




Figure 5 Binding of 125 l-labeled CryICa to S. exigua BBMV. (A) Homologous competition. 125 l-labeled Cryl C was incubated with BBMV from 
susceptible (FRA) or resistant (Xen-R) insects in the presence of increasing concentrations of unlabeled CrylC toxin. (B) Dissection of total binding 
into their reversible and irreversible components. Contribution of both components to the binding of 125 l-labeled CryICa to BBMV from 
susceptible and resistant insects were compared. T-test {P = 0.7244 for the reversible component) {P = 0.0170 for the irreversible component). 
BBMV, brush border membrane vesicles. 
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Figure 6 Effect of suppression of SeABCC2 or SeABCC3 gene expression on pupation (A, B, and C) and adult emergence (D) of 

S. exigua. The expression of SeABCC2 or SeABCC3 gene was suppressed by feeding dsSeABCC2 or dsSeABCC3 to newly molted second instar 

larvae. A viral gene (ORF302)-specific dsRNA (dsCON) was used for a dsRNA control. Each treatment was replicated three times. (A) Pupation, 

(B) pupal weight, (C) pupal duration and (D) adult emergence after dsRNA feeding were measured. Each replicate used 10 larvae. Different letters 

above standard deviation bars indicate significant differences among means at Type I error = 0.05 (LSD test). dsRNA, double-stranded RNA, LSD, 

least significant difference. 



discriminate susceptible insects in the BC offsprings we 
have only used one of the bulk extremes (the resistant in- 
sects). As a consequence, evidence of linkage to resistance 
was not deduced by the comparison of the frequency of 
the markers in the two bulks, but by comparing their 
frequency deviation from the expected 50% and their 
deduced genomic localization. 2) Due to the absence of 
genomic information on S. exigua, the resistance gene has 
been positioned and deduced using the genome of a re- 
lated Lepidoptera, B. mori. 3) Genetic markers have been 
discovered by RNA sequencing instead of direct examin- 
ation of genomic DNA. BSA based on RNA information 
could introduce a certain bias because the cDNA contri- 
bution to the bulk can be allele specific. To prevent or 
minimize such a possibility the signal has been detected at 
the genomic level and not at the gene level, thus averaging 
out the possible bias present in one particular gene. Our 



approach using NGS and BSA applied to S. exigua has 
pinpointed two genomic regions with high segregation 
distortion that have a high probability of containing genes 
involved in the resistance to B. thuringiensis. BSA has 
been previously employed for the identification of a pesti- 
cide resistance gene [30]; however, as far as we know, this 
is the first report of successful mapping of insect resist- 
ance genes by combining asymmetric BSA with RNA 
sequencing. 

One of these two regions contained three genes cod- 
ing for ABC type C transporters, including the ABCC2 
gene along with ABCC1 and ABCC3 genes. The se- 
quences of the ABCC2 and ABCC3 genes coded for 12- 
transmembrane domain proteins. Most interestingly, 
the ABCC2 allele from the resistant insects contained a 
deletion that eliminates a cytoplasmic portion of the 
protein affecting one of the two ATP binding regions, 
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making this gene a candidate for the resistance in the 
Xen-R colony. The second region, located in chromo- 
some 22, did not reveal genes potentially involved in the 
mode of action of B. thuringiensis, although we cannot 
discard the presence of additional genes in such a re- 
gion contributing to the resistance in the Xen-R insects. 
Alternatively, synteny is not complete in Lepidoptera 
and it could be possible that both regions are located in 
the same chromosome in S. exigua. 

To determine the role of the ABCC2 transporter, and 
that of its closest paralog ABCC3, in the mode of action of 
Cryl proteins, their developmental expression and tissue 
distribution has been studied. They both were expressed 
in the larval midgut, in agreement with the expected char- 
acteristics from membrane proteins that must interact 
with Cryl proteins toxic to larvae. However, at early stages 
of development, when insects are more susceptible to B. 
thuringiensis, the expression levels of ABCC3 are lower 
than those of ABCC2, suggesting that the latter may play 
a more relevant role in the mode of action of B. thurin- 
giensis toxins. A more direct approach to the involvement 
of ABCC2 and ABCC3 in the toxicity of Cry proteins was 
to measure the susceptibility to Cryl Ac and CrylCa in lar- 
vae where the expression of these two transporters had 
been decreased by feeding dsRNA. With both toxins, mor- 
tality significantly decreased when the expression of either 
ABCC2 or ABCC3 was reduced. This result shows that 
both transporters are involved in the toxicity of these two 
Cryl proteins and that, in theory, mutations in either one 
could result in resistance to these two toxins. Since the se- 
quence homology between ABCC2 and ABCC3 is higher 
than 70%, it is likely that both proteins could serve as re- 
ceptors for Cryl toxins. ABCC1, although not as similar 
to ABCC2 as ABCC3, cannot be discarded at this point as 
also playing a role in the mode of action of B. thuringien- 
sis. However, the most direct evidence of the involvement 
of ABCC transporters in the mode of action of B. thur- 
ingiensis toxins in S. exigua is that resistant insects, but 
not susceptible controls, carried a mutated allele of the 
ABCC2 gene. 

All previous reports relating the ABCC2 transporter to 
resistance to CrylA proteins have found mutations in the 
ABCC2 gene that affected the extracellular domains in dif- 
ferent ways. In H. virescens, the mutation consisted of a 
22-bp deletion in exon 2, which would result in a trun- 
cated product of 99 amino acids instead of the 1,339 
amino acids in the wild type protein [20]. In P. xylostella, 
the ABCC2 gene had a 30-bp deletion in exon 20, which 
was predicted to remove the final transmembrane domain 
with the result of positioning the C-terminal end of the 
protein in the extracellular space instead of its natural lo- 
cation inside the cell [21]. In B. mori, the mutation added 
an extra tyrosine residue in the second extracellular loop 
of the ABCC2 transporter [22]. So far, the mutation in the 



ABCC2 gene conferring resistance in T. ni has not yet 
been elucidated [21]. Binding studies with the above resist- 
ant colonies showed that insects with the ABCC2 muta- 
tion lacked binding of CrylAb and CrylAc in the case of 
H. virescens and P. xylostella [20,31]. The same lack of 
binding of these two toxins was found in the T. ni resistant 
colony [32]. However, no apparent alteration of CrylAb 
binding (CrylAc was not tested) was found in the B. mori 
resistant insects as compared to susceptible controls [22]. 
As has been hypothesized elsewhere [24], it is very likely 
that the extra tyrosine residue in the extracellular loop 
does not prevent binding of CrylAb, although it could 
confer resistance by altering the transporter function and 
decreasing the probability of the toxin being drawn into 
the membrane. 

The mutation we have found in the S. exigua resistant 
insects is the first reported mutation that does not affect, 
directly or indirectly, the extracellular ABCC2 domains, 
but affects the intracellular C-terminal part of the pro- 
tein, which contains the second ATP binding domain. If 
the above model of the involvement of the ABCC2 
transporter in the mode of action of Cryl toxins is cor- 
rect, we should not expect any decrease in toxin binding 
since the extracellular domains of the protein are not al- 
tered, and the intracellular deletion does not seem to ap- 
preciably affect the conformation of the rest of the 
protein. Similarly to B. mori, resistance would be con- 
ferred by altering the transporter function affecting its 
role in inserting the toxin into the membrane [24]. Our 
binding data completely agree with this hypothesis. First, 
no significant difference was found regarding the affinity 
of the toxin to bind BBMV, and only a small difference 
was found between resistant and susceptible insects re- 
garding the concentration of binding sites. Second, alter- 
ation of the transporter ability to draw the toxin into the 
membrane would specifically affect the irreversible com- 
ponent of binding, and this is precisely what our results 
indicate. The involvement of the ABCC2 transporter in 
the resistance to B. thuringiensis in S. exigua adds a new 
species in which the role of this gene in resistance to 
Cry proteins has been demonstrated. 

Silencing of ABCC2 and ABCC3, in addition to allow- 
ing us to confirm the role of these receptors in the 
mode of action of B. thuringiensis, has enabled the study 
of the fitness cost associated with the absence of these 
proteins. Although dsRNA feeding could not achieve 
complete knock-out of gene expression to emulate a re- 
sistance gene, partial silencing revealed clear effects on 
insect fitness. Fitness cost has been associated with re- 
sistance in most of the B. thuringiensis resistant strains 
[33]. According to our results, field selection of ABC- 
related resistance alleles can be counteracted by the 
reduced fitness of resistant insects versus susceptible in- 
sects in areas not exposed to B. thuringiensis products 
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(that is, refuges) and, in this way, contribute to delay the 
development of resistance. 

Previous studies on the Xen-R strain revealed that the 
insects, in the absence of exposure to B. thuringiensis, had 
a high level of expression of certain genes {arylphorin and 
repatS among them) that are up-regulated in susceptible 
insects after B. thuringiensis intoxication [25]. For this 
reason we found it interesting to determine whether a de- 
crease in the expression, or a malfunction, of the SeABCC2 
or SeABBC3 transporters could trigger the up-regulation 
of such genes. Our results reveal that silencing of 
SeABCC2 or SeABCC3, in susceptible insects, induces the 
expression of arylphorin and repatS. Such a result sug- 
gests that higher expression levels of these genes in the 
Xen-R insects are induced by the mutation in the 
SeABCC2 receptor. Additionally, up-regulation of these 
genes in response to B. thuringiensis intoxication might be 
induced by the interaction of the toxins with the ABCC 
receptors, perhaps because binding of the toxin would 
probably block their physiological function. Other studies 
have shown that arylphorin and repatS, as well as other 
members of the repat superfamily, are induced in response 
to intoxication with several types of toxins in Spodoptera 
spp. [34-36]. It would be of interest to check whether 
other resistant strains with mutations in the ABCC2 trans- 
porter overexpress these genes and to explore the possibil- 
ity of using such a phenotype as a marker to monitor for 
the presence of ABC-related mutations in the field. 

Also worth noting is the fact that the ABCC2 mutation 
in S. exigua apparently confers resistance not only to 
CrylA proteins but to CrylCa as well. Although Xentari™ 
contains CrylA and CrylCa as major Cry components 
[37], S. exigua larvae are susceptible to CrylCa but rela- 
tively tolerant to CrylA proteins [38]. Therefore, resist- 
ance to Xentari™ in Xen-R implies resistance to the 
CrylCa key component. Another indirect piece of evi- 
dence that mutations in the ABCC2 transporter may be 
involved in resistance to Cryl toxins other than those of 
the CrylA family was recently reported in P. xylostella 
[39]. CrylA resistant insects with a mutation in the 
ABCC2 gene were cross-resistant to Cryl Fa and they 
lacked binding of CrylA proteins and Cryl Fa to BBMV 
[39]. In this insect species, it was already known that 
Cryl Fa shares binding sites with CrylA proteins [40]. 

Fifteen ABCC genes have been annotated on the 
B. mori genome [41,42], It would be interesting to 
test whether additional members of such a family in 
Lepidoptera, but also in other insect orders, could be 
involved in modulating the insect susceptibility to Cryl 
toxins and also to other groups of toxins, such as Cry2 
and Cry3. Examples of membrane proteins that are 
known to serve as surrogate receptors for a wide variety 
of Cry proteins are cadherin, APN and ALP. All have 
been reported to be the binding proteins for quite 



different Cry toxins, such as Cryl, Cry3, and Cry4 and 
Cryll, which are active against Lepidoptera, Coleoptera 
and Diptera, respectively [10,43-48]. It seems that dif- 
ferent Cry proteins have found the appropriate insect 
species, or even insect order, to interact with a same 
type of surrogate receptor. 

Conclusions 

BSA, based on high-throughput transcriptome sequen- 
cing, can be applied for genetic mapping in inver- 
tebrates. In this study, BSA has contributed to the 
identification of a novel type of mutation in the ABCC2 
transporter conferring high levels of insect resistance to 
B. thuringiensis. Additional results have expanded the 
role of the ABCC2 transporter in B. thuringiensis resist- 
ance beyond the CrylA family of proteins to include 
CrylCa and revealed the involvement of different mem- 
bers of the ABCC family in the mode of action of B. 
thuringiensis and in the fitness-cost associated with the 
resistance to B. thuringiensis. 

Methods 

Source of Cry toxins 

B. thuringiensis Cry toxins, Cryl Ac and CrylCa, were ob- 
tained from recombinant Escherichia coli XL-1 strains 
(kindly supplied by Dr. Ruud de Maagd, Plant Research 
International, Wageningen, The Netherlands). Protein pro- 
duction, solubilization and trypsin activation was performed 
as described earlier [49]. Proteins used for bioassays were in 
the form of protoxin dissolved in carbonate buffer (0.1 M 
Na 2 C0 3 , 0.01 M ethylenediaminetetraacetic acid (EDTA), 
0.17 M NaCl, pH 10.5). For binding assays, CrylCa was fur- 
ther purified by anion-exchange chromatography followed 
by gel filtration on a Superdex-75 column containing 
1 mM dithiothreitol (DTT) according to published proce- 
dures [50,51], to remove the adsorbed toxin fragments that 
could interfere in the radiolabeling of the toxin. Protein con- 
centration was measured using the Coomassie protein assay 
reagent (Bio-Rad, Hercules, CA, USA) with bovine serum al- 
bumin (Sigma- Aldrich, Tres Cantos, Spain) as the standard. 

Source of insects and genetic crosses 

For genetic and binding analyses, two S. exigua laboratory 
colonies were used. The resistant colony (Xen-R) was 
highly resistant to the commercial bioinsecticide Xentari™ 
(based on B. thuringiensis subsp. aizawai, Valent Biosci- 
ences, Libertyville, IL, USA). This colony has been previ- 
ously described by Hernandez-Martinez et al, [25] and the 
resistance level was maintained with a constant selection 
protocol in which neonate larvae were reared for five days 
on an artificial diet containing Xentari™ [52], at 10 mg per 
gram of artificial diet. The susceptible control colony 
(FRA) originated from France and has been maintained 
for at least 10 years without B. thuringiensis exposure. For 
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the functional analysis of ABCC genes, a third S. exigua 
colony (KOR) was established from larvae collected from a 
field population infesting welsh onion in Andong, Korea. 
Larvae from all colonies were reared on an artificial diet at 
25°C and adults were fed with 10% sucrose. Insects were 
reared on an artificial diet at 25 ± 3°C with 70 ± 5% relative 
humidity and a photoperiod of 16/8 hours (light/dark). 

For resistance mapping, single-pair crosses were per- 
formed between females from the Xen-R colony and 
males from the FRA colony. From the Fl offspring from 
each family, one male was backcrossed with five Xen-R 
females. Neonate larvae from the resulting progeny (BC 
families) were reared for five days on an artificial diet 
containing Xentari™ (10 mg/g) and then transferred to a 
normal diet until they reached the fifth instar. Survivors 
from each BC family were pooled (Sel-BC sample) and 
used for NGS and BSA. Sel-BC insects must be homozy- 
gous for the resistance allele or alleles and it should con- 
tain, on average, a 75% genome background of the 
parental resistant colony. 

RNA extraction and NGS sequencing 

For BSA, three samples were processed: Xen-R, FRA, 
and Sel-BC. For each sample, about 30 early fifth instar 
larvae midguts were dissected and used for total RNA 
isolation using Trizol (Invitrogen, Carlsbad, CA, USA) 
following the manufacturers instructions. The quality 
and quantity of total RNA was determined spectroscop- 
ically at 230, 260, and 280 nm. About 20 \ig of total 
RNA from each sample was sent to Macrogen Inc 
(Seoul, Korea) for mRNA isolation and high-throughput 
sequencing using the Illumina Hi-Seq2000 paired-end 
technology. 

Sequence analysis 

Hi-Seq2000 reads were analyzed using seq_scrumbs [53] 
to trim bad quality regions. These reads were assembled 
with Trinity [54] along with the 454-sequencing reads of 
the previous transcriptome [28]. To keep only the most 
informative unigenes, a selection of the assembled uni- 
genes was carried out. For a unigene to be finally 
selected, it should have to be evaluated as positive for 
any one of the following criteria: having a CDS (Coding 
DNA sequence), having a Blastx match against Bombyx 
or Drosophila proteins, being expressed with a level 
higher than 0.25 RPKM (reads per kilo base per million) 
or being longer than 500 bp. 

A SNVs (single nucleotide variants: SNPs and Indels) 
search was done only on the longest unigene from each 
defined unigene cluster. The Illumina and 454-sequencing 
reads were mapped against these reduced unigene sets 
with bowtie2 [55] and bwa for sequences from the 454- 
sequencing [56]. SNVs were identified with Varscan2 soft- 
ware [57]. Orthologues for each unigene were identified 



with a bi-directional Blast search comparison against B. 
mori and Drosophila melanogaster protein databases. 

Analysis of the genomic region 

The genomic sequence for the ABCC2 region, which in- 
cludes the mutation in the Xen-R colony, was obtained. 
Primers flanking the deleted region were designed based 
on the mRNA sequence from the FRA colony and used 
for PCR amplification of the genomic region using total 
DNA from Xen-R and FRA adults. PCR products for 
both colonies were cloned in the pGEM-T-easy vector 
(Invitrogen) and several clones were sequenced and as- 
sembled to generate the consensus sequence. At least 
three independent clones were used for each of the se- 
quenced samples. 

Sequence comparison and phylogenetic analysis 

Sequences from the ABCC receptors from other insect spe- 
cies used for phylogenetic analyses were retrieved after 
Blast-search against different databases: NCBI, Silkdb and 
Manduca Base (Agripest). For the sake of clarity, only rep- 
resentative genes from selected lepidopteran species and 
from the coleopteran T. castaneum were included in the 
analysis. Amino acid sequences were aligned using ClustalX 
with default settings [58] and visualized in GeneDoc [59]. 
Phylogenetic analyses were performed using the neighbor- 
joining method with 100 bootstraps, using the ClustalX 
[58] and MEGA5 programs [60]. Phylogenetic trees were 
visualized with the MEGA5 program. Nucleotide binding 
domains in the ABCC receptors were identified by com- 
parison against the CDD (conserved domain database) at 
NCBI [61]. Transmembrane regions in the protein were 
predicted using TOPCONS online software [62]. 

Semiquantitative and quantitative RT-PCR 

Total RNA was extracted from whole bodies of different 
developmental stages or from different tissues from fifth 
instar larvae with Trizol (Invitrogen) according to the 
manufacturers instructions. Aliquots (1 \ig) were incu- 
bated at 70°C for three minutes and then used for syn- 
thesizing cDNA using an RT-mix kit (Intron, Seoul, 
Korea). For semiquantitative RT-PCR, the synthesized 
single-stranded cDNA was used as a template for PCR 
amplification with gene-specific primers for SeABCC2 
(5 -GAG CCT AAC ATG GCT GAC TAT C-3' and 5 - 
CTG CCT ACA TTT GCG TCT ACT-3 ) and SeABCC3 
(5' -TCT ATA CGA CCG CTC TCT TAC C-3' and 5 - 
CAC CCT GCA CAC CTT CTT AT-3 ) with 35 cycles 
under the following conditions: 20 seconds at 94°C for 
denaturation, 50 seconds at 55°C for annealing, and one 
minute at 72°C for extension. A housekeeping gene, 
actin, was used as an internal control in each sample for 
an equivalent of template with actin primers (5'-TGG 
CAC CAC ACC TTCTAC-3' and 5 -CAT GAT CTG 
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GGT CAT CTT CT-3) with 35 cycles under the follow- 
ing conditions: 20 seconds at 94°C for denaturation, 
20 seconds at 50°C for annealing, and one minute at 72°C 
for extension. 

The mRNA expression level of SeABCC2 and SeABCC3 
genes was further confirmed by quantitative RT-PCR 
(qRT-PCR). qRT-PCR was performed with an Applied Bio- 
systems 7500 real-time PCR system (Applied Biosystems, 
Foster City, CA, USA) using SYBR® Green Realtime PCR 
master mix (Toyobo, Osaka, Japan) according to the 
manufacturers instructions with the gene-specific primers 
described above. Template cDNA from 48 hour dsRNA- 
treated third instar larvae (see below) was used for the 
qRT-PCR reaction. After a hot start at 94°C for 10 minutes, 
qRT-PCR was performed with 40 cycles of one minute at 
94°C, 30 seconds at 55°C and 40 seconds at 72°C. The 
housekeeping gene fi-actin was used as a control with the 
primers described above. Each condition was independ- 
ently replicated three times. Quantitative analysis of 
SeABCC2 and SeABCC3 expression was done using the 
comparative C T (AAC T ) method [63]. 

Preparation and delivery of RNA interference (RNAi) of 
SeABCC2 and SeABCC3 

Template cDNA was amplified with the following primers 
containing T7 RNA polymerase promoter sequence for 
SeABCC2 (5 -TAA TAC GAC TCA CTA TAG GGA 
GAG GAG CCT AAC ATG GCT GAC TAT C-3' and 5 - 
TAA TAC GAC TCA CTA TAG GGA GAG CTG CCT 
ACA TTT GCG TCT ACT-3 ) and SeABCC3 (5 -TAA 
TAC GAC TCA CTA TAG GGA GAG TCT ATA CGA 
CCG CTC TCT TAC C-3' and 5 -TAA TAC GAC TCA 
CTA TAG GGA GAG CAC CCT GCA CAC CTT CTT 
AT-3'). Sequence identity between SeABCC2 and SeABCC3 
at the region used for dsRNA was less than 65%. The PCR 
products were used to prepare dsRNA using the MEGA 
Script RNAi kit according to the manufacturers instruc- 
tions (Ambion, Austin, TX, USA). The synthesized RNAs 
were annealed at 37°C for four hours and then left at 70°C 
for five minutes. Before dsRNA feeding, the solution was 
prepared by mixing a transfection reagent, Metafectene 
PRO (Biontex, Plannegg, Germany), with dsRNA in 1:1 
(v/v) ratio and then incubated for 20 minutes at 25°C. 

For the dsRNA feeding assays, cabbage (lxl cm) was 
used as a carrier by applying 4 |iL of the dsRNA (150 ng 
total) solution on a piece of cabbage leaf. Second instar 
larvae were starved for six hours before being allowed to 
feed on cabbage treated with dsRNA for four hours at 
25°C. Then, fresh artificial diet was supplied to the lar- 
vae to continue development. To determine the effect of 
the transfection reagent on S. exigua larvae, a control 
bioassay was performed using only 4 \iL of Metafectene 
PRO solution on cabbage leaf. Metafectene PRO does 
not affect the development and mortality of control 



larvae (data not shown). Knockdown of SeABCC2 and 
SeABCC3 gene expression was evaluated by RT-PCR at 
selected periods up to 96 hours post feeding. A viral 
gene, CpBV-ORF302, was used as a negative control for 
RNAi [64]. 

Bioassays 

S. exigua third instar larvae were used in bioassays with 
CrylAc and CrylCa proteins. For the concentration re- 
sponse bioassays, the concentrations tested for CrylAc 
were 0.625, 1.25, 2.5, 5 and 10 [ig/cm 2 , while for CrylCa, 
the concentrations tested were 0.062, 0.125, 0.25, 0.5 and 
1 (ig/cm 2 . Cryl toxins were serially diluted with sterile de- 
ionized water and then overlaid on the piece of cabbage 
leaf and air-dried for 20 minutes. Each bioassay was con- 
ducted with 10 larvae per replicate and three replicates 
per concentration. Mortality was scored at five days post- 
treatment. LC 50 values were calculated by Probit analysis 
using the EPA Probit Analysis program, version 1.5 (US 
Environmental Protection Agency, Cincinnati, OH, USA). 

For dsRNA bioassays, S. exigua second instar larvae 
were fed dsRNA specific to SeABCC2 and SeABCC3 as 
described above for 24 hours. Then, larvae were exposed 
to CrylAc or CrylCa pro toxins for five days using the 
same method as before. Control larvae were treated with 
4 \iL of control dsRNA (dsCON). Each treatment was 
replicated three times with ten larvae per replicate. 
Means and variances of treatments were analyzed in 
one-way analysis of variance (ANOVA) by PROC GLM 
of SAS program (SAS Institute, 1989). The means were 
compared by the LSD test at Type I error = 0.05. 

Binding assays 

Last-instar larvae from the susceptible (FRA) and the re- 
sistant (Xen-R) colonies were dissected and the bolus- 
free midguts were preserved at -80°C until required. 
BBMVs were prepared by the differential magnesium 
precipitation method [65], frozen in liquid nitrogen, and 
stored at -80°C. The total protein concentration of the 
BBMV preparations was determined by the Bradford 
method using bovine serum albumin as standard. 

Chromatography-purified CrylCa was labeled using 
the chloramine-T method [66]. CrylCa (25 \ig) was 
mixed with 0.5 mCi of 125 NaI (PerkinElmer, Boston, 
MA, USA), chloramine T was added to be at 6 mM final 
concentration. After 45 seconds of incubation at room 
temperature, the reaction was stopped by adding potas- 
sium metabisulfite followed by non- radioactive Nal. The 
specific activity of the 125 I-CrylCa was 0.73 \iCU\ig. 

All binding assays were performed at room temperature 
in a final volume of 0.1 ml in binding buffer (phosphate 
buffer saline: pH 7.4 phosphate-buffered saline (PBS), sup- 
plemented with 0.1% bovine serum albumin). For the de- 
termination of the optimal concentration of BBMV to be 
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used for the binding assays, 1 nM of I-CrylCa was in- 
cubated for one hour with increasing concentrations of 
BBMV (from either FRA or Xen-R insects). The non- 
specific binding was determined using an excess of un- 
labeled toxin. 

Competition binding assays were performed by incu- 
bating 1 nM of 125 I-CrylCa with 32 [ig of BBMV and in- 
creasing concentrations of unlabeled toxin for one hour. 
The unbound radio-labeled toxin was washed out and 
the amount of residual bound 125 I-CrylCa was deter- 
mined using a 2480 WIZARD 2 (Perkin Elmer, Waltham, 
MA, USA) gamma counter. The dissociation constant 
(K d ) and concentration of binding sites (R t ) were calcu- 
lated using the LIG AND program [67] . 

To determine the contribution of reversible and irre- 
versible binding to the observed specific binding, three re- 
action mixtures were prepared with 1 nM of 125 I-CrylCa 
and 30 \ig of BBMV, and all the samples were incubated 
for two hours. One sample was used to determine the total 
binding. In a second sample, used to determine the non- 
specific binding, an excess of unlabeled toxin was added at 
the start of the incubation. To the third sample, which 
was used to measure the irreversible binding, an excess of 
unlabeled toxin was added after one hour of incubation 
and the reaction was allowed to proceed until two hours. 
The specific binding was calculated by subtracting the 
radioactivity in the pellet of the second sample from that 
in sample one and the irreversible binding was calculated 
by subtracting the radioactivity in the pellet of the second 
sample from that in the third sample. Therefore, the 
reversible binding was calculated by subtracting the irre- 
versible binding from the amount of specific binding. Ex- 
periments were performed in triplicate and similar results 
were obtained from two independent batches of BBMV. 

Developmental analysis of dsRNA treatment 

RNAi-treated second instar larvae were analyzed for the 
possible effects on pupation and adult emergence. Meta- 
fectene PRO transfection reagent (4 \iL) without dsRNA 
was used to feed control larvae. Each treatment was repli- 
cated three times with ten larvae per replicate. Pupation 
success, weight of newly pupated individuals, pupation 
time at 25°C and adult emergence were analyzed. 

Expression levels of arylphorin and repat5 genes by RNAi of 
SeABCC2 and SeABCC3 genes 

cDNAs prepared from individuals at 48 hours after 
dsSeABCC2, dsSeABCC3 or dsRNA 0 *™ 2 (dsCON) feed- 
ing were used for the analysis of arylphorin and repatS 
gene expression by qRT-PCR. After a hot start at 94°C 
for 10 minuts, qRT-PCR amplification was performed as 
before with gene-specific primers for arylphorin and 
repatS [25]. 
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